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ABSTRACT 

In  this  study,  the  propagation  and  growth  of  whistler  mode  waves 
generated  by  electron  beams  within  Garth's  bow  shock  is  investigated 
using  a  planar  model  for  the  bow  shock  and  a  model  electron 
distribution  function.  Within  the  shock,  the  model  electron 
distribution  function  possesses  a  field-aligned  T^  >  T#  beam  that  is 
directed  toward  the  magnetosheath.  Waves  with  frequencies  between 
about  1  and  100  Hz  with  a  wide  range  of  wave  normal  angles  are 
generated  by  the  beam  via  Landau  and  anomalous  cyclotron  resonances. 
However,  because  the  growth  rate  is  small  and  because  the  wave  packets 
traverse  the  shock  quickly,  these  waves  do  not  attain  large  amplitudes. 
Waves  with  frequencies  between  about  30  and  150  Hz  with  a  wide  range  of 
wave  normal  angles  are  generated  by  the  beam  via  the  normal  cyclotron 
resonance.  The  ray  paths  for  most  of  these  waves  are  directed  toward 
the  solar  wind  although  some  wave  packets,  due  to  plasma  convection, 
travel  transverse  to  the  shock  normal.  These  wave  packets  grow  to 
large  amplitudes  because  they  spend  a  long  time  in  the  growth  region. 
The  results  suggest  that  whistler  mode  noise  within  the  shock  should 
increase  in  amplitude  with  increasing  upstream  3gn.  The  study  provides 
an  explanation  for  the  origin  of  much  of  the  whistler  mode  turbulence 
observed  at  the  bow  shock. 
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I.  INTRODUCTION 

The  discovery  of  field-aligned  electron  beams  within  the  Earth's 
bow  shock  by  the  ISEE-2  fast  plasma  experiment  [Feldman  et  al.,  1982] 
is  an  important  contribution  to  existing  measurements  of  electron  dis¬ 
tribution  functions  at  collisionless  shocks.  Two  important  features  of 
the  beams  are  that  they  are  moving  toward  the  magnetosheath  and  that 
the  beam  energy,  which  ranges  up  to  about  100  eV,  increases  with 
increasing  distance  from  the  solar  wind  boundary  of  the  shock. 

Although  the  observed  characteristics  of  the  electron  beams  have  not 
yet  been  categorized  with  respect  to  upstream  solar  wind  parameters,  it 
is  known  that  the  beams  are  observed  for  both  quasi-parallel  and  quasi¬ 
perpendicular  bow  shock  geometries.  It  is  thought  that  the  beams  are 
produced  by  a  macroscopic  electric  field  within  the  shock.  To  obtain  a 
qualitative  understanding  of  how  this  electric  field  may  be  created, 
the  reader  can  consult  the  paper  by  Wu  [1982]. 

It  is  natural  that  the  discovery  of  field-aiigned  electron  beams 
at  the  bow  shock  would  initiate  studies  on  instabilities  generated  by 
the  beams.  The  interest  in  plasma  instabilities  stems  from  the  fact 
that  many  types  of  plasma  waves  are  often  observed  at  the  bow  shock 
(e.g.,  Fredricks  et  al.  [1968]  and  Rodriguez  and  Gurnett  [1975]).  In 
addition,  plasma  turbulence  is  thought  to  play  an  important  dissipative 
role  at  the  shock.  Since  the  discovery  of  the  beams,  two  studies  have 
investigated  the  linear  Vlasov  stability  of  the  observed  electron 


4 


distributions.  In  Thomsen  et  al.  (1983]  it  is  shown  that  the  electron 
beams  generate  electrostatic  waves  via  Landau  resonance  while  in  Tokar 
et  al.  [1984]  it  is  shown  that  the  beams  can  also  generate  whistler 
mode  waves  via  Landau  and  cyclotron  resonance.  Because  this  study  is 
an  extension  of  the  work  reported  in  Tokar  et  al.  [1984],  it  is  desir¬ 
able  to  briefly  discuss  the  main  results  of  that  paper. 

Figure  1,  reproduced  from  Tokar  et  al.  [1984],  illustrates  typical 
observations  of  whistler  mode  waves  at  the  bow  shock.  In  the  top  panel 
observations  of  the  wave  magnetic  fields,  measured  by  the  ISEE  1  plasma 
wave  experiment  on  December  13,  1977,  are  shown  in  ten  frequency 
channels  covering  the  range  5.6  Hz  to  1  kHz.  The  bottom  panel  depicts 
the  magnetic  field  magnitude  measured  by  the  magnetic  field  experiment 
on  ISEE  1.  The  magnetic  ramp  defining  the  shock  transition  region  is 
clearly  visible,  with  the  upstream  boundary  at  about  17:35:05.0. 

Within  the  magnetic  ramp,  whistler  mode  noise  can  be  seen  in  the  bottom 
six  frequency  channels  of  the  wave  magnetic  field  data.  The  noise  is 
most  intense  within  the  shock  and  extends  into  the  upstream  solar  wind 
and  downstream  magnetosheath. 

To  investigate  the  possibility  that  the  field-aligned  electron 
beams  present  In  the  shock  generate  the  whistler  mode  noise,  Tokar  et 
al.  [1984]  analyzed  the  linear  Vlasov  stability  of  three  analytic  fits 
to  electron  velocity  distribution  data  obtained  by  the  fast  plasma 
experiment  within  the  shock.  The  analysis  shows  that  the  three  elec¬ 
tron  distribution  functions  are  unstable  to  whistler  mode  waves.  In 
all  cases,  the  free  energy  regions  in  the  distribution  function  that 


are  responsible  for  the  wave  growth  are  attributable  to  the  field- 
aligned  beam.  Waves  are  generated  by  Landau  and  cyclotron  resonance 
and  have  wave  vectors  directed  toward  both  the  magnetosheath  and  the 
solar  wind.  The  unstable  waves  have  plasma  rest  frame  frequencies 
between  about  1  and  100  Hz  and  wave  normal  angles  that  range  from  0°  to 
the  resonance  cone  angle.  In  addition,  the  growth  rates  for  many  of 
the  waves  appeared  to  be  large  enough  to  give  rise  to  the  large  ampli¬ 
tude  whistler  mode  noise  shown  in  Figure  1. 

This  study  extends  the  analysis  given  in  Tokar  et  al.  [1984]  and 
has  two  primary  purposes.  The  first  is  to  investigate  the  generation 
of  whistler  mode  waves  throughout  the  shock  transition  region,  and  the 
second  is  to  investigate  the  propagation  and  path  integrated  amplifica¬ 
tion  of  the  waves.  The  previous  study  treated  only  three  positions 
within  the  shock.  In  this  study,  a  simple  model  electron  distribution 
function  is  constructed  that  is  a  function  of  position  within  the  shock 
transition  region.  The  distribution  evolves  from  an  isotropic 
Maxwellian  at  the  solar  wind  to  an  isotropic  modified  Lorentzian  at  the 
magnetosheath  and  generates  a  broad  spectrum  of  whistler  mode  waves. 
Using  a  simple  model  for  the  shock,  ray  paths  for  the  generated  whis¬ 
tler  mode  waves  are  calculated.  The  path  integrated  amplification 
along  the  ray  path  is  obtained  using  the  model  electron  distribution 
function  and  a  linear  growth  rate  formula. 


II.  MODEL  BOW  SHOCK  AND  ELECTRON  DISTRIBUTION  FUNCTION 


In  this  section  the  simple  models  for  the  bow  shock  region  and  the 
electron  distribution  function  are  presented.  These  models  provide  the 
description  of  the  geometry  and  plasma  parameters  within  the  shock  that 
is  needed  to  investigate  the  generation  and  propagation  of  the  whistler 
mode  waves. 


A,  Planar  Bow  Shock  Model 

For  simplicity,  a  planar  model  is  adopted  to  describe  the  bow 
shock  region.  The  planar  assumption  is  justified  because  the  length 
scales  along  the  shock  normal  and  transverse  to  the  shock  normal  are 
both  much  less  than  the  radius  of  curvature  of  the  shock.  Figure  2 
illustrates  the  coordinate  system  and  overall  geometry  of  the  model. 
The  z  axis  is  along  the  shock  normal  with  positive  z  directed  from  the 
solar  wind  to  the  magnetosheath.  The  shock  transition  region  extends 
for  100  km  in  the  z  direction  and  separates  the  upstream  solar  wind 
from  the  downstream  magnetosheath.  The  chosen  thickness  of  the  shock, 
100  km,  is  in  good  agreement  with  the  measurements  of  Russell  et  al. 

[ 1979J . 

In  order  to  calculate  ray  paths  for  whistler  mode  waves  using  the 
planar  bow  shock  model,  the  electron  concentration,  magnetic  field 
vector,  and  plasma  flow  vector  must  be  known  throughout  the  shock 
transition  region.  In  the  planar  model,  these  quantities  are  a 


function  only  of  z,  the  penetration  into  the  shock.  They  are  assumed  to 
vary  as  follows: 

1)  The  electron  concentration,  n,  and  magnetic  field  magnitude, 

B,  increase  linearly  within  the  shock  transition  region  from 
upstream  solar  wind  values  to  downstream  magnetosheath  values. 
The  increase  in  n  across  the  shock  is  from  2.5  cm-3  to  10  cm-3 
while  the  increase  in  B  across  the  shock  is  from  5  gammas  to 
20  gammas.  These  values  correspond  to  the  values  measured 
during  the  December  13,  1977,  shock  crossing  analyzed  in  Tokar 
et  al.  [1984]. 

2)  The  direction  of  the  magnetic  field  follows  from  the  assumed 
increase  in  the  magnitude  of  if  and  the  condition  that  the 
component  of  if  in  the  direction  of  the  shock  normal,  n,  remain 
constant  throughout  the  shock  transition  region.  Figure  2 
illustrates  a  magnetic  field  line  for  the  case  when  the 
upstream  angle  between  if  and  n,  0an,  is  45°.  Two  values  of  0gn 
are  treated,  45°  and  65°. 

■+  +> 

3)  The  plasma  flow  velocity,  Vf,  is  taken  to  be  in  the  z  direction 
throughout  the  shock  transition  region.  The  magnitude  of  vj 
decreases  with  increasing  z  such  that  the  normal  component  of 
the  number  flux,  n|vf|,  remains  constant  throughout  the  transi¬ 
tion  region.  In  the  model  adopted,  the  plasma  flow  speed 
decreases  from  400  km/s  at  the  solar  wind  to  100  km/s  at  the 
magnetosheath. 

In  summary,  within  the  shock  transition  region  the  electron  con¬ 
centration  and  magnetic  field  magnitude  increase  linearly,  the  normal 


component  of  B  is  constant,  the  plasma  flow  velocity  is  in  the  z 
direction,  and  the  normal  component  of  the  number  flux  is  constant. 

It  should  be  noted  that  the  adopted  increases  in  the  electron 
concentration  and  magnetic  field  magnitude  across  the  shock,  intended 
to  agree  with  observation,  do  not  satisfy  MHD  jump  conditions  for 
shocks.  A  model  based  on  the  MHD  jump  conditions  can  be  constructed, 
for  example,  by  reducing  the  jump  in  the  electron  concentration  until  a 
desired  jump  in  the  magnetic  field  magnitude  is  obtained.  The  cal¬ 
culated  upstream  and  downstream  parameters  must  then  be  connected  in 
some  fashion  across  the  shock.  The  differences  in  the  adopted  model 
and  a  model  constructed  using  the  MHD  jump  conditions  do  not  signifi¬ 
cantly  effect  the  primary  results  of  this  study.  In  particular,  the 
neglect  of  a  transverse  plasma  flow  within  the  shock  does  not 
significantly  affect  the  integrated  growth  calculations  because  the 
wave  packets  that  grow  to  large  amplitudes  have  vg  primarily  In  the 
transverse  direction  and  |vg|  >  |vf|. 

It  is  also  evident  that  familiar  features  at  supercritical  shocks, 
such  as  overshoots  in  |iT|  and  p  and  a  magnetic  "foot",  are  not  included 
in  the  model  adopted.  Such  features  can  affect  both  the  calculated  ray 
paths  and  growth  rates.  However,  it  is  difficult  to  construct  a  model 
electron  distribution  function  throughout  the  shock  that  reflects  these 
structures.  For  this  reason,  the  simple  planar  model  with  linearly 
increasing  |£j  and  n  is  adopted  to  obtain  a  semi -quantitative  picture 
of  the  propagation  of  the  whistler  mode  waves. 


B.  Model  Electron  Distribution  Function 


It  is  desirable  to  obtain  a  model  electron  distribution  function 
that  is  a  function  of  position  within  the  shock  transition  region  so 
that  the  generation,  propagation,  and  growth  (or  damping)  of  whistler 
mode  waves  throughout  the  transition  region  can  be  studied  in  detail. 

In  this  study,  a  simple  model  electron  distribution  function  is  con¬ 
structed  such  that  its  features  are  consistent  with  experimental 
results.  An  alternative  approach  would  be  to  solve  the  Vlasov-Maxwell 
equations  throughout  the  shock.  While  such  a  solution  would  possess  a 
self-consistency  absent  in  the  simple  model  adopted,  obtaining  an 
accurate  solution  would  be  a  difficult  task.  Consequently,  the  program 
adopted  is  to  construct  a  simple  model  that  agrees  with  experimental 
findings. 

The  analytic  function  that  often  provides  a  good  fit  to  the  elec¬ 
tron  distribution  functions  measured  in  the  shock  is  the  sum  of  a 
modified  Lorentzian  and  a  convected  Maxwellian  [Feldman  et  al.,  1983]. 
Specifically,  the  electron  distribution  function  data  in  the  plasma 
rest  frame  is  fit  to  a  function  of  the  form 

F(vi*v«)  ”  Fj  +  F2  (1) 


where 
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is  the  modified  Lorentzian  and 
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F2  *  C2exp ^  2kT^  )  exp^2kT^'*  (3) 

is  the  convected  Maxwellian.  In  these  equations,  and  vj  are  the 
electron  speeds  perpendicular  and  parallel  to  the  magnetic  field, 
respectively,  and  T n  are  the  perpendicular  and  parallel  temperatures 
of  the  beam,  respectively,  vj  is  the  breakpoint  speed  of  the  modified 
Lorentzian,  V2  is  the  convection  speed  of  the  beam,  m  is  the  electron 
mass,  k  is  Boltzmann's  constant,  and  C},  C2 ,  R  and  P  are  constants. 

Referring  to  Equations  2  and  3  it  can  be  seen  that  the  following 
eight  parameters  must  be  known  throughout  the  shock:  Cj ,  C2,  vj,  V2, 

T II ,  T j_,  P  and  R.  To  construct  functions  that  describe  the  variation  of 
these  parameters  with  position  in  the  shock  transition  region  the 
following  experimental  results  are  employed. 

1)  As  the  transition  region  Is  traversed  by  the  spacecraft,  the 
electron  distribution  function  evolves  smoothly  from  the 
upstream  solar  wind  distribution  to  the  downstream  magneto¬ 
sheath  distribution.  A  typical  solar  wind  electron  distribu¬ 
tion  function  is  an  isotropic  Maxwellian  with  a  temperature  of 
about  1.5  *  105°K.  A  typical  magnetosheath  electron  distribu¬ 
tion  function  is  a  modified  Lorentzian  with  R  *  4.0,  P  *  R-l 
and  V[  «  6  *  10®  cm/s  (Feldman  et  al.,  1983]. 


2)  The  electron  distribution  function  in  the  shock  transition 


region  has  both  a  beam  component  and  a  flat-topped  component. 
The  convection  speed  of  the  beam,  V2»  increases  with 
increasing  shock  penetration.  The  increase  is  consistent  with 
the  interpretation  that  the  incoming  solar  wind  electrons  are 
accelerated  by  a  macroscopic  electric  field  within  the  shock 
transition  region.  The  breakpoint  speed  of  the  flat-topped 
component  of  the  distribution  is  approximately  equal  to  the 
beam  convection  speed  (i.e.,  vj  -  V£)  [Feldman  et  al.,  1983]. 

3)  The  perpendicular  temperature  of  the  beam,  T^,  is  accurately 
modeled  by  assuming  that  the  first  adiabatic  invariant  of  the 
electrons  is  conserved  as  the  electrons  move  from  the  solar 
wind  to  the  magnetosheath.  The  parallel  temperature  of  the 
beam,  Tj,  decreases  with  increasing  shock  penetration.  The 
decrease  in  Tj  is  qualitatively  explained  by  the  fact  that 
the  electrons  are  being  accelerated  along  ?.  A  value  of  T g 
equal  to  about  2  x  10l+oK  is  typical  near  the  center  of  the 
shock  transition  region  [Feldman  et  al.,  1983]. 

Consistent  with  the  above  points,  the  following  simple  models  are 
chosen  to  describe  the  evolution  of  the  electron  distribution  function 
within  the  shock  transition  region. 

In  agreement  with  point  1  above,  the  distribution  evolves  from  an 
isotropic  Maxwellian  with  a  temperature  of  1.5  x  105°K  at  the  solar 
wind  (z  M  0  km)  to  an  isotropic  modified  Lorentzian  with  vj  ■  6  *  108 
cm/s  and  R  ■  4,  P  ■  R  -  1  at  the  magnetosheath  (z  ■  100  km).  This 


evolution  is  accomplished  through  the  coefficients  Cj  and  C2  in  Equa¬ 
tions  2  and  3.  The  coefficient  Cj  in  Equation  2  increases  linearly 
from  the  value  0.  at  the  solar  wind  to  1.  at  the  magnetosheath  while  C2 
in  Equation  3  decreases  linearly  from  the  value  1.  at  the  solar  wind  to 
0.  at  the  magnetosheath.  An  additional  constant  multiplies  the  entire 
electron  distribution  function  to  ensure  that  the  integral  of  F(v^,vj) 
over  all  velocity  space  is  equal  to  the  assumed  electron  concentration. 
This  constant  is  a  function  of  z  because,  in  the  planar  bow  shock 
model,  the  electron  concentration  is  a  function  of  z.  It  can  be 
verified  that  at  about  z  ■  30  km  half  of  the  particles  reside  in  the 
convected  Maxwellian  and  half  reside  in  the  modified  Lorentzian. 

In  agreement  with  point  2  above,  vj  is  set  equal  to  V2»  The  con¬ 
vection  speed  of  the  beam,  V2»  ranges  from  0  cm/s  at  z  -  0  km  to  about 
6  x  108  cm/s  at  z  =  100  km.  As  a  function  of  z,  the  convection  speed 
increases  as  would  the  speed  of  an  electron  that  is  being  accelerated 
by  a  constant  electric  field  of  magnitude  1  volt/km.  This  corresponds 
to  a  potential  drop  across  the  shock  of  100  volts  and  results  in  beam 
energies  within  the  shock  transition  region  that  are  in  good  agreement 
with  the  measurements  of  Feldman  et  al.  [1983], 

In  agreement  with  point  3  above ,  the  perpendicular  temperature  of 
the  beam  is  calculated  assuming  that  the  first  adiabatic  invariant  of 
the  electrons  is  conserved  throughout  the  transition  region.  The 
parallel  temperature  of  the  beam  decreases  linearly  from  1.5  x  105oK  at 
z  -  0  km  to  2  x  104oK  at  z  -  50  km.  For  z  greater  than  50  km,  T y  is 


equal  to  2  x  104°K 


The  above  discussion  defines  the  electron  distribution  function 
throughout  the  shock  transition  region.  Contour  plots  of  the  base  10 
logarithm  of  the  model  distribution,  spanning  one  order  of  magnitude  of 
F(vj.,V||),  are  shown  in  Figure  3  for  six  positions  within  the  shock 
transition  region.  This  figure  illustrates  that  the  distribution  func¬ 
tion  evolves  from  a  relatively  cool  Maxwellian  at  the  solar  wind  (z  *  0 
km)  to  a  hotter  modified  Lorentzian  at  the  magnetosheath  (z  ■  100  km). 
With  increasing  z,  the  convection  speed  and  ratio  Tj/Tj  for  the  beam 
increase. 


III.  WHISTLER  MODE  STABILITY  OF  THE  MODEL  DISTRIBUTION 


Combining  the  model  electron  distribution  function  described  in 
the  previous  section  with  an  expression  for  the  linear  growth  rate  of 
whistler  mode  waves,  the  growth  of  whistler  mode  waves  throughout  the 
shock  transition  region  can  be  investigated.  The  equation  used  to 
calculate  the  growth  rate  is  identical  to  the  equation  given  in  Kennel 
and  Wong  [1967]  and  Tokar  et  al.  [1984].  Because  a  detailed  discussion 
of  the  whistler  mode  instability  mechanisms  for  electron  distributon 
functions  measured  in  the  shock  transition  region  is  given  in  Tokar  et 
al.  [1984],  only  a  brief  discussion  will  be  given  here.  The  primary 
purpose  of  this  section  is  to  summarize  the  unstable  wave  frequencies 
for  the  model  electron  distribution  function. 

For  the  model  electron  distribution  function  described  in  the 
previous  section,  the  unstable  whistler  mode  waves  frequencies  and  wave 
normal  angles  are  a  function  of  penetration  into  the  transition  region, 
z.  This  is  because  the  evolution  with  z  of  the  velocity  distribution 
causes  the  conditions  for  instability  to  be  functions  of  z.  Before 
presenting  the  unstable  frequencies  for  the  model  distribution,  a  brief 
discussion  of  how  the  evolution  of  the  model  distribution  with  z 
affects  the  conditions  for  instability  is  desirable. 

In  Tokar  et  al.  [1984]  it  Is  shown  that  for  generation  via  the 
normal  and  anomalous  cyclotron  resonances  the  distribution  must  possess 
regions  of  sufficiently  large  anisotropy  at  the  cyclotron  resonance 


velocities.  In  the  model  distribution  this  anisotropy  is  due  to  the 
increase  in  the  beam  temperature  ratio,  T x/T  N »  with  increasing  z. 
Because  the  cyclotron  resonance  velocities  must  remain  in  regions  of 
large  anisotropy  and  because  the  beam  energy  increases  with  z,  as 
illustrated  in  Figure  3,  it  is  expected  that  the  unstable  frequencies 
will  be  functions  of  z.  In  Tokar  et  al.  [1984]  it  is  also  shown  that 

for  growth  via  the  Landau  resonance,  the  distribution  must  possess  a 

k  || 

positive  and  sufficiently  large  -r-i — r  -r- —  at  the  wave  phase  velocity. 

|k,||  dv  u 

From  an  inspection  of  Figure  3  it  is  clear  that  as  the  Maxwellian 
component  of  the  model  distribution  evolves  away  from  the  vj  **  0  axis, 
this  condition  will  be  satisfied  for  k H  <  0  if  z  is  less  than  about  50 
km.  As  z  increases  from  50  km  to  100  km,  the  parallel  derivative  of 
F(vi,V|)  decreases  to  zero  because  the  dominant  component  of  the  dis¬ 
tribution  is  the  isotropic  modified  Lorentzian. 

Figures  4  and  5  summarize  the  unstable  plasma  rest  frame  frequen¬ 
cies  for  the  model  electron  distribution  function  described  in  the  pre¬ 
vious  section.  It  should  be  noted  that  the  unstable  frequencies  are 
shown  only  as  a  function  of  position  within  the  shock.  The  unstable 
wave  normal  angles  for  a  given  frequency  are  not  summarized.  This  is 
because  the  results  obtained  in  Tokar  et  al.  [1984]  apply  to  the  model 
electron  distribution  function  with  the  main  point  being  that  for  a 
given  frequency  an  appreciable  range  of  wave  normal  angles  is  usually 
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In  Figure  4  the  unstable  frequencies  for  generation  via  the  normal 
cyclotron  resonance  are  shown  as  a  function  of  z.  This  resonance  is 
the  m  *  -1  contribution  to  the  growth  rate  expression  employed  by  Tokar 
et  al.  [1984].  For  this  generation  mechanism,  the  waves  all  have 
kq  >  0,  corresponding  to  wave  vectors  directed  toward  the  solar  wind. 
The  spectrum  of  unstable  frequencies  is  in  good  agreement  with  the 
range  obtained  for  the  three  electron  distributions  analyzed  in  Tokar 
et  al.  [1984],  The  increase  in  the  magnitudes  of  the  unstable  frequen¬ 
cies  with  increasing  z  locates  the  normal  cyclotron  resonance  velocity 
in  the  region  of  large  anisotropy  associated  with  the  field-aligned 
beam.  The  distribution  is  stable  for  small  values  of  z  because  the 
beam  temperature  ratio,  Tj/T|,  is  small.  Near  the  center  of  the  shock 
transition  region  the  cyclotron  resonant  interactions  are  strongest. 
This  is  because  in  this  region  Tj/T(|  is  large  and  the  Maxwellian  compo¬ 
nent  of  the  model  distribution  is  appreciable. 

In  Figure  5  the  unstable  frequencies  for  generation  via  the  Landau 
and  anomalous  cyclotron  resonances  are  shown  as  a  function  of  z.  These 
resonances  are  the  m  ■  0  and  +1  contributions  to  the  growth  rate 
expression  employed  by  Tokar  et  al.  [1984],  For  these  generation 
mechanisms,  the  waves  all  have  kj  <  0,  corresponding  to  wave  vectors 
directed  toward  the  magnetosheath.  In  agreement  with  expectations,  the 
distribution  is  stable  for  values  of  z  greater  than  about  50  km.  The 
range  of  unstable  frequencies  is  in  good  agreement  with  the  range 
obtained  for  the  distributions  analyzed  in  Tokar  et  al.  [1984]. 
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IV.  RAY  PATHS  AND  INTEGRATED  AMPLIFICATION 


In  the  last  section,  the  unstable  whistler  mode  frequencies  are 
summarized  for  the  model  electron  distribution  function  as  a  function 
of  position  within  the  shock  transition  region.  In  this  section  the 
subsequent  propagation  and  growth  (or  damping)  of  the  generated 
whistler  mode  waves  are  discussed. 

Figures  4  and  5  show  that  within  the  shock  transition  region  the 
growth  rate  for  a  given  whistler  mode  wave  frequency  is  positive  in 
only  a  limited  region.  The  positions  within  the  shock  where  the  growth 
rate  is  positive  act  as  amplification  regions  for  whistler  mode  wave 
packets.  To  compute  the  overall  amplification  we  assume  that  a  back¬ 
ground  of  very  weak  whistler  mode  noise  with  a  wide  range  of  wave 
normal  angles  is  available  for  amplification  from  the  solar  wind. 
Neubauer  et  al.  [1977]  have  reported  that  a  background  continuum  of 
whistler  mode  noise  is  a  nearly  ubiquitous  feature  of  the  solar  wind. 

In  this  study  the  ray  paths  emanating  from  the  source  location  are 
calculated  using  the  planar  bow  shock  model  and  the  WKB  approximation 
(i.e.,  the  geometrical  optics  limit).  The  amplification  along  the  ray 
path  is  calculated  using  the  model  electron  distribution  function  and 
the  linear  growth  rate  expression  given  in  Kennel  and  Wong  [1967]  and 
employed  in  Tokar  et  al.  [1984],  Although  the  validity  of  this  proce¬ 
dure  in  the  turbulent  shock  transition  region  may  be  questionable  at 
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low  frequencies  where  Che  wavelength  becomes  comparable  Co  Che  thick¬ 
ness  of  Che  transition  region,  Che  results  provide  at  least  a  qualita¬ 
tive  understanding  of  the  subsequent  propagation  of  the  generated 
waves . 

A  first-order  validity  criterion  for  ray  tracing  in  the  transition 
region  is  that  the  characteristic  scale  length  for  changes  in  the 
plasma  parameters  must  be  larger  than  the  wavelength  of  the  waves 
investigated.  A  more  quantitative  validity  criterion  for  a  magnetized 
plasma  is  difficult  to  derive  [Budden,  1961],  Because  the  character¬ 
istic  scale  length  is  the  thickness  of  the  transition  region,  100  km, 
in  this  study  the  ray  tracing  calculations  are  limited  to  waves  with 
wavelengths  less  than  about  25  km.  This  condition  eliminates  consider¬ 
ation  of  the  low  frequency  (1-20  Hz)  waves  generated  by  the  Landau  and 
anomalous  cyclotron  resonances.  However,  this  restriction  is  not 
serious  because  the  results  indicate  that  only  the  waves  generated  by 
the  normal  cyclotron  resonance  are  amplified  to  significant  levels 
above  the  background  noise. 

To  facilitate  the  understanding  of  the  ray  tracing  procedure,  the 
steps  followed  in  a  typical  calculation  will  be  listed  individually. 
Before  proceeding,  it  should  be  noted  that  the  procedure  is  complicated 
somewhat  by  the  fact  that  the  plasma  flow  speed  is  a  function  of  the 
position  coordinate  within  the  shock,  z.  The  situation  is  analogous  to 
the  refraction  of  sound  waves  produced  by  a  gradient  of  the  wind  velo¬ 
city  (e.g.,  see  Landau  and  Lifshitz  [1959]).  Although  this  effect  is 
included  in  the  calculations,  the  change  in  the  wave  normal  angle  along 
the  ray  path  due  to  the  gradient  in  the  plasma  flow  speed  is  small 


compared  to  the  changes  produced  by  the  gradients  in  the  electron 
concentration  and  magnetic  field  magnitude. 


The  steps  followed  in  the  ray  tracing  procedure  are  given 
below. 


1)  For  a  chosen  position  within  the  shock  transition  region,  z, 
a  real  wave  frequency,  ujj.,  and  initial  wave  normal  angle,  0O, 
are  selected  for  which  the  growth  rate,  u^(z),  is  positive. 

In  this  study  the  initial  wave  normal  angles  0°  and  45°  are 
analyzed.  For  the  initial  wave  normal  angle  of  0°,  the  direc¬ 
tion  of  the  wave  vector  is  along  the  magnetic  field.  However, 
for  the  Initial  wave  normal  angle  of  45°,  the  possible  wave 
vectors  lie  on  a  cone  with  axis  the  magnetic  field  direction 
and  half  angle  equal  to  45°.  Of  this  set  of  initial  wave 
vector  orientations  the  two  wave  vectors  that  lie  in  the  plane 
of  the  planar  bow  shock  model  are  investigated. 

2)  The  magnitude  of  the  wave  vector,  |ic|  ,  and  the  group  velocity 
of  the  wave  packet ,  Vg ,  are  calculated  using  the  wave  normal 
angle  and  the  cold  plasma  whistler  mode  dispersion  relation. 
For  the  frequencies  of  interest  in  this  study,  a  good  approxi¬ 
mation  for  the  dispersion  relation  is  [Stix,  1964] 
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In  this  equation,  u>p  is  the  electron  plasma  frequency,  is 
the  electron  gyrof requency ,  and  c  is  the  speed  of  light.  The 
values  of  cop  and  o)g  at  position  z  are  calculated  using  the 
magnetic  field  magnitude  and  electron  concentration  models 
adopted  for  the  planar  shock  model. 

3)  To  obtain  the  ray  path  in  the  frame  of  reference  of  the  planar 
bow  shock  model,  the  wave  packet  is  then  advanced  along  the 
vector  Vg  +  vf(z)  for  a  fixed  increment  of  z,  Az.  The  vector 
Vf(z)  is  the  plasma  flow  vector  at  the  position  z  in  the 
planar  bow  shock  model. 

4)  The  amplification  along  the  segment  of  the  ray  path  connecting 
z  and  z  +  Az  is  then  calculated.  Denoting  the  length  of  the 
segment  by  As,  the  amplification  along  the  segment,  Ar,  is 
approximated  by 

0)  (z)  As 

Ar(z)-_i - - -  .  (5) 

|vg(z)  +  vf(z) | 

Over  this  segment  of  the  ray  path,  the  fields  are  amplified  by 
the  factor  e^. 

5)  Next,  the  wave  normal  angle  9(z  +  Az)  is  calculated  at  the 
position  z  +  Az.  To  do  so,  first  write  the  Doppler  relation 
for  o)'(z  +  Az),  the  wave  frequency  at  z  +  Az  in  the  frame  of 
reference  of  the  plasma  parcel  at  z.  This  equation  is 
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where  Av  -  Vf(z  +  Az)  -  Vf(z)  is  the  relative  velocity  of  the 
plasma  parcels  located  at  z  and  z  +  Az.  Next,  combine  Equa¬ 
tions  4  and  6  to  obtain  the  magnitude  of  the  wave  vector  at 
z  +  Az 


|k(z+Az)  | 
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a)  cos  9(z+Az)  -  a)(z)  -  k(z+Az)«Av 
8 


(7) 


In  the  planar  bow  shock  model,  Snell's  law  requires  that  the 
component  of  ic  perpendicular  to  the  shock  normal  remain 
constant  along  the  ray  path.  Because  it  is  difficult  to  solve 
Equation  7  for  |k(z  +  Az) j ,  in  the  right-hand  side  of  Equation 
7  the  approximation  is  made  that  tc(z)  »  li(z  +  Az).  This 
approximation  retains  the  £•  Av  correction  to  the  wave  fre¬ 
quency  and  is  accurate  for  small  increments  in  z.  With  this 
approximation,  Equation  7  is  solved  for  |&(z  +  Az) | .  Snell's 
law  then  yields  the  equation  needed  to  calculate  the  new  wave 
normal  angle,  9(z  +  Az). 

6)  Finally,  steps  2  through  5  are  repeated  with  the  new  fre¬ 
quency,  o)'(z  +  Az),  and  wave  normal  angle,  9(z  +  Az).  The 
procedure  is  continued  until  the  net  amplification  along  the 
ray  path,  T,  which  is  the  sum  of  the  values  of  Ar  for  each  ray 
path  segment,  becomes  negative. 


Figures  6  and  7  illustrate  the  relative  orientations  of  the  mag¬ 
netic  field  vector,  the  wave  vector,  the  group  velocity,  and  the  plasma 
flow  velocity  within  the  shock  transition  region.  Figure  6  applies  to 
a  wave  packet  that  is  generated  by  the  normal  cyclotron  resonance 
(m  =  -1  and  ky  >  0)  while  Figure  7  applies  to  a  wave  packet  that  is 
generated  by  the  Landau  and  anomalous  cyclotron  resonances  (m  ■  0,  +1 
and  k||  <  0).  The  two  initial  wave  vector  orientations  for  the  wave 
normal  angle  0  are  shown.  The  orientations  are  distinguished  by  the 
angle  with  the  <p  =  0°  orientation  corresponding  to  the  wave  vector 
with  the  largest  component  in  the  direction  of  the  solar  wind.  The  <{>  =* 
180°  orientation  lies  on  the  cone  of  half  angle  0  and  is  in  the  plane 
of  the  planar  bow  shock  model.  This  orientation  is  denoted  <|>  *  180° 
because  it  can  be  visualized  as  being  180°  in  azimuth  on  the  cone  from 
the  4>  =  0°  orientation. 

The  group  velocity  Vg,  the  plasma  flow  velocity,  Vf,  and  the  angle 
between  Vg  and  B,  are  also  shown  in  Figures  7  and  8.  The  relative 
lengths  of  vg  and  vf  shown  in  these  figures  are  typical  for  calcula¬ 
tions  near  the  center  of  the  shock  transition  region.  Figure  6 
illustrates  that  when  k n  is  positive  the  plasma  convection  causes  the 
ray  paths  to  be  refracted  transverse  to  the  shock  normal  while  Figure  7 
illustrates  that  when  ky  is  negative  the  plasma  convection  causes  the 
ray  paths  to  be  refracted  predominantly  in  the  direction  of  the  magnet¬ 
osheath. 

The  results  of  the  ray  tracing  calculations  are  summarized  in 
Figures  8  and  9.  Figure  8  applies  to  the  case  when  the  upstream  angle 
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between  B  and  n,  0]jn,  in  the  planar  bow  shock  model  is  equal  to  45° 
while  Figure  9  is  for  0jjn  equal  to  65°.  These  two  values  of  0gn 
illustrate  the  important  effects  that  occur  as  the  shock  becomes  more 
perpendicular.  Parallel  or  nearly  parallel  geometries  and  perpendi¬ 
cular  or  nearly  perpendicular  geometries  are  not  treated  in  this  study 
due  to  a  lack  of  information  concerning  the  electron  distribution 
function  for  these  shocks. 

The  magnetic  field  line  and  several  ray  paths  within  the  shock 
transition  region  are  shown  in  Figure  8.  The  ray  paths  are  labelled 
using  wave  frequency,  initial  wave  normal  angle,  and  value  of  <J>  in  the 
format  (u>r/2n,  0O,  $).  A  value  of  <p  for  0  *  0°  is  not  applicable  and 
so  is  omitted  from  the  ray  path  labels.  The  solid  portions  of  the  ray 
paths  correspond  to  regions  of  growth  while  the  dashed  portions  corre¬ 
spond  to  regions  of  damping.  The  amplification  factor,  T,  at  the  posi¬ 
tion  of  marginal  stability,  *  0,  is  given  at  the  heads  of  the  arrows 
denoting  the  ray  path  directions. 

The  ray  paths  for  generation  via  the  Landau  and  anomalous  cyclo¬ 
tron  resonances  lie  near  the  center  of  the  transition  region  in  Figure 
8.  The  frequency  analyzed,  u^./2tt  =  50  Hz,  has  associated  wavelengths 
always  less  than  25  km  and  serves  to  illustrate  the  important  propaga¬ 
tion  characteristics  for  the  case  k y  <  0.  Two  initial  wave  vector 
orientations,  0O  ■  45°  with  <p  ■  0°  and  0O  =  45°  with  <p  -  180°  are 


treated.  The  case  0O  ■  0°  is  omitted  from  the  calculations  because  the 
growth  rate  is  always  negative  for  k j  <  0  and  0O  *  0°  [Tokar  et  al., 
1984). 
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Figure  8  shows  that  the  ray  paths  for  the  frequency  0^/2 tt  *■  50  Hz 
are  directed  toward  the  magnetosheath.  This  is  due  to  the  fact  that 
both  the  group  velocity  and  the  solar  wind  velocity  are  directed  toward 
the  magnetosheath,  as  shown  in  Figure  7.  For  generation  via  the  Landau 
and  anomalous  cyclotron  resonances,  the  net  amplification  along  the  ray 
paths  is  small  because  the  waves  propagate  quickly  through  the  transi¬ 
tion  region.  The  values  of  T  given  in  Figure  8,  0.5  and  0.9,  are  typi¬ 
cal  for  all  wave  packets  with  k.  u  <  0. 

The  ray  paths  for  generation  via  the  normal  cyclotron  resonance 
are  shown  originating  on  the  right-hand  side  of  Figure  8.  The  three 
frequencies  treated,  0^/2^  =  70,  90,  and  110  Hz,  together  with  the 
three  initial  wave  vector  orientations,  illustrate  all  of  the  important 
propagation  characteristics  for  the  case  k#  >  0.  Figure  7  illustrates 
that  in  this  case  the  ray  paths  are  directed  toward  the  solar  wind  and 
the  amplification  factors  can  be  large.  In  none  of  the  cases  do  the 
wave  packets  reach  the  solar  wind. 

Figure  9  illustrates  the  results  of  the  ray  tracing  calculations 
when  0j}n  in  the  planar  bow  shock  model  is  equal  to  65°.  The  format  of 
the  figure  is  identical  to  that  in  Figure  8.  The  new  magnetic  field 
geometry  does  not  significantly  affect  the  calculations  for  the  waves 
traveling  toward  the  magnetosheath.  Again,  the  primary  result  is  that 
the  background  noise  is  not  significantly  amplified. 

However,  the  new  magnetic  field  geometry  does  affect  the  calcula¬ 
tions  for  the  waves  generated  by  the  normal  cyclotron  resonance.  As  is 
clear  from  an  inspection  of  Figure  9,  many  of  the  ray  paths  for  genera¬ 
tion  via  the  normal  cyclotron  resonance  have  large  amplification 


factors.  In  particular,  the  largest  amplification  occurs  when  the  ray 
path  is  directed  nearly  transverse  to  the  shock  normal.  These  waves 
grow  to  large  amplitudes  because  they  spend  a  long  time  in  the  growth 
region.  In  two  cases,  labeled  (70,  45,  180)  and  (110,  45,  180),  the 
waves  asymptotically  approach  a  point  where  the  plasma  flow  velocity 
and  the  z  component  of  the  group  velocity  are  of  equal  magnitude  and 
oppositely  directed,  in  a  region  where  the  growth  rate  is  positive. 
Consequently,  these  wave  packets  do  not  leave  the  growth  region  and 
grow  to  Infinite  amplitude.  This  infinite  growth  is  denoted  in  the 
figure  by  the  symbol  °°.  For  the  same  two  frequencies  but  different 
source  regions,  the  amplification  factors  did  not  return  to  zero  within 
1000  km  of  the  source  region.  These  ray  paths  are  terminated  at  the 
point  of  marginal  stability,  and  the  values  for  T  at  the  point  of 
marginal  stability,  13.3  and  140.7,  are  given. 

It  should  be  noted  that  if  components  of  it  and  Vg  out  of  the  plane 
of  the  bow  shock  model  are  considered,  the  additional  degree  of  freedom 
will  result  in  additional  wave  packets  that  satisfy  Vg »z  equal  to 
-|vf|.  Therefore,  a  more  accurate  three-dimensional  model  will  possess 
more  wave  packets  attaining  large  amplitude,  thereby  supporting  the 
conclusions  of  this  study. 

The  results  of  the  ray  tracing  calculations  are  summarized  as 
follows. 

1)  The  waves  generated  by  the  Landau  and  anomalous  cyclotron 

resonances  follow  ray  paths  that  are  directed  toward  the  mag¬ 
netosheath.  The  wave  packets  traverse  the  transition  region 


quickly,  experiencing  a  period  of  small  growth  followed  by  a 
period  of  small  damping.  Because  these  waves  propagate 
quickly  through  the  transition  region  they  are  not  signifi¬ 
cantly  amplified.  Large  amplitude  wave  packets  with  k|  <  0 
are  therefore  not  expected  to  reach  the  magnetosheath. 

2)  Many  of  the  waves  generated  by  the  normal  cyclotron  resonance 
achieve  large  amplitudes.  The  waves  with  largest  amplitude 
lie  in  the  frequency  range  from  about  70  to  110  Hz  and  are 
generated  near  the  center  of  the  shock  transition  region,  in  a 
region  where  the  normal  component  of  the  propagation  velocity 
(vf  +  Vg)  is  small. 

3)  As  the  shock  classification  changes  from  quasi-parallel  to 
quasi-perpendicular,  the  intensity  of  whistler  mode  noise  at 
the  shock  increases.  This  is  because  the  time  a  wave  packet 
spends  in  the  growth  region  increases  with  0gn.  This  effect 
has  been  observed  at  the  bow  shock  [Rodriguez  and  Gurnett, 


V.  DISCUSSION 


The  integrated  growth  calculations  presented  in  the  previous  sec¬ 
tion  suggest  that  the  electron  beams  present  in  the  bow  shock  generate 
whistler  mode  waves  with  large  amplification  factors  and  that  the  waves 
are  absorbed  within  the  shock  transition  region.  At  this  point  it  is 
desirable  to  compare  the  predictions  of  the  study  with  observations  of 
whistler  mode  noise  at  the  shock.  Two  comparisons  can  be  made.  The 
predicted  frequencies  of  the  waves  can  be  compared  with  observation  and 
the  measured  amplitudes  of  the  waves  can  be  compared  with  the  inte¬ 
grated  growth  calculations.  However,  because  the  initial  angular  dis¬ 
tributions  of  the  observed  waves  are  not  known,  the  predicted  wave 
vector  directions  cannot  be  compared  with  observation. 

As  was  discussed  in  the  last  section,  the  instabilities  generated 
by  the  normal  cyclotron  resonance  have  the  largest  integrated  growth 
rates.  For  this  reason,  the  unstable  frequencies  for  generation  via 
this  resonance  should  be  compared  to  the  frequencies  for  the  large 
amplitude  noise  shown  in  Figure  1.  In  Figure  1,  the  range  in  frequency 
over  which  the  whistler  mode  noise  is  observed  is  about  5.6  Hz  to  100.0 
Hz.  From  the  previous  section,  the  predicted  plasma  rest  frame  fre¬ 
quencies  for  the  large  amplitude  waves  generated  via  the  normal  cyclo¬ 
tron  resonance  range  from  about  70  Hz  to  110  Hz.  This  leads  to  the 
conclusion  that  the  two  frequency  ranges  are  in  rough  agreement. 


This  agreement  is  improved  if  the  Doppler-shif t  caused  by  the 
relative  motion  of  the  satellite  and  the  plasma  is  taken  into  account. 
Because  the  wave  vectors  of  many  of  the  waves  generated  by  the  normal 
cyclotron  resonance  are  directed  into  the  solar  wind  flow,  the  Doppler- 
effect  for  these  waves  results  in  a  measured  frequency  that  is  lower 
than  the  plasma  rest  frame  frequency.  When  the  wave  vector  has  a 
component  in  the  direction  of  the  solar  wind  flow,  the  measured  fre¬ 
quency  will  be  higher  than  the  plasma  rest  frame  frequency.  It  is  not 
necessary  to  work  out  the  Doppler-shif t  for  all  rays  because  a  typical 
calculation  will  illustrate  the  magnitude  of  the  effect  involved. 

The  Doppler-shifted  frequency  of  the  wave  can  be  calculated  using 
the  equation 

in  *  oj'  +  £  »v^  (8) 


where  ui  is  the  Doppler-shifted  frequency,  u>'  is  the  rest  frame  fre¬ 
quency,  lc  is  the  rest  frame  wave  vector,  and  Vf  is  the  plasma  flow 
speed.  To  illustrate  the  general  effect,  consider  a  rest  frame  fre- 
quency  100  Hz,  the  two  wave  vector  orientations  of  9  ■  45°  with  <p  m  0° 
and  0  ■  45°  with  <J>  *  180°,  and  the  case  9gn  *  45°.  From  the  results 
presented  in  Figure  8,  these  wave  packets  exist  near  z  -  50  km  and  have 
large  integrated  growth  rates.  Employing  Equation  8  to  calculate  the 
measured  frequency  for  the  two  wave  vector  orientations,  it  is  found 
that  for  the  <p  *  0°  orientation  the  frequency  is  downshifted  to  about 


1 


91  Hz  and  for  Che  $  ■  180°  orientation  the  frequency  is  upshifted  to  to 
about  105  Hz. 

This  Doppler  shift  can  increase  dramatically  if  highly  oblique 
wave  normal  angles  are  considered.  This  is  because  the  magnitude  of 
the  wave  vector,  |ic|,  becomes  large  as  the  wave  normal  angle  approaches 
the  resonance  cone  angle,  %g8  =  cos-^^/^).  However,  because  hot 
plasma  effects  on  the  whistler  mode  dispersion  relation  may  be  impor¬ 
tant  near  the  resonance  cone,  this  possibility  will  not  be  investigated 
in  detail  here. 

It  is  also  important  to  investigate  whether  or  not  the  instabilit¬ 
ies  of  this  study  can  attain  the  large  amplitudes  of  the  shock  asso¬ 
ciated  whistler  mode  noise  shown  in  Figure  1.  From  Neubauer  et  al. 
[1977],  it  is  found  that  typical  spectral  densities  of  whistler  mode 
noise  In  the  solar  wind  near  1  AU  are  10“3  to  10-2  gammas/ (Hz)  1  /2  for 
the  frequency  7  Hz  and  10”1*  gammas/ (Hz)  1/2  for  the  frequency  70  Hz. 
These  amplitudes  for  the  background  noise  are  to  be  compared  to  ampli¬ 
tudes  of  whistler  mode  noise  in  the  shock  transition  region.  An  esti¬ 
mate  of  the  amplification  required  can  be  obtained  from  the  ISEE-1 
plasma  wave  receiver  on  December  13,  1977.  Figure  1  illustrates  that 
at  17:35:15.0  UT  on  this  day  ISEE-1  is  in  the  shock  transition  region 
and  large  amplitude  whistler  mode  noise  is  present.  At  this  time, 
spectral  densities  for  the  noise  are  about  1  gamma/(Hz)  1/2  at  7  Hz  and 
10"2  gamma/ (Hz) 1 /2  at  70  Hz. 


From  these  results,  it  Is  concluded  that  the  total  amplification 
of  the  background  noise  must  be  a  factor  of  102  or  10^.  This  corre¬ 
sponds  to  a  range  in  T  from  about  5  to  7.  It  is  clear  from  an  inspec¬ 
tion  of  Figures  8  and  9  that  many  of  the  waves  generated  via  the  normal 
cyclotron  resonance  have  amplification  factors  large  enough  to  give 
rise  to  the  large  amplitude  whistler  mode  noise  at  the  shock.  This 
noise  is  expected  to  have  maximum  intensity  near  the  center  of  the 
shock  transition  region  or  magnetic  ramp.  This  is  in  agreement  with 
Figure  1  where  the  noise  has  maximum  intensity  near  the  center  of  the 
shock  transition  region  at  about  17:35.13.0  UT. 

This  study  demonstrates  that  many  of  the  whistler  mode  waves 
generated  by  electron  beams  in  the  bow  shock  can  attain  the  large 
amplitudes  typical  of  observations  at  the  shock.  The  calculated  fre¬ 
quencies  of  the  waves  and  the  observed  frequencies  of  whistler  mode 
noise  at  the  shock  are  in  rough  agreement.  The  agreement  is  Improved 
if  the  Doppler-shifts  of  oblique  whistlers  are  taken  into  account.  The 
results  of  the  study  also  suggest  that  the  amplitudes  of  the  whistler 
mode  noise  at  the  shock  should  Increase  with  increasing  6gn,  in  agree¬ 
ment  with  the  observations  of  Rodriguez  and  Gurnett  [1975J. 


The  results  of  the  study  do  not  provide  a  direct  explanation  for 
the  origin  of  the  whistler  mode  noise  in  the  solar  wind  and  magneto¬ 
sheath  boundary  regions  of  the  shock.  It  is  plausible  that  this 
discrepancy  is  a  consequence  of  the  idealized  models  adopted  and  that 
in  actuality  most  of  the  whistler  mode  noise  observed  at  the  shock  is 
generated  in  the  transition  region  and  propagates  into  the  upstream 
solar  wind  and  downstream  magnetosheath.  For  example,  if  the  wave 
packets  depicted  in  Figures  8  and  9  encounter  electron  distributions 
near  the  solar  wind  or  magnetosheath  that  are  weakly  or  marginally 
stable  to  whistler  mode  waves,  they  will  escape  the  transition  region. 
This  idea  is  consistent  with  a  large  number  of  plasma  wave  observations 
at  shocks.  In  many  cases,  whistler  mode  turbulence  is  most  intense  at 
the  shock  and  decreases  in  intensity  approximately  exponentially  with 
increasing  distance  from  the  shock. 
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FIGURE  CAPTIONS 


Figure  1  This  figure  is  reproduced  from  Tokar  et  al.  [1984] 

and  illustrates  typical  observations  of  whistler  mode 
noise  at  the  bow  shock.  The  figure  depicts  data 
obtained  with  the  ISEE-1  plasma  wave  experiment  and 
the  magnetic  field  experiment.  The  top  panel  shows 
the  wave  magnetic  fields  in  ten  frequency  channels 
while  the  bottom  panel  illustrates  the  magnetometer 
data.  The  magnetic  ramp  defining  the  shock  transi¬ 
tion  region  is  clearly  visible  as  is  the  large  ampli¬ 
tude  whistler  mode  noise  at  the  shock.  The  noise  is 
most  intense  within  the  shock  transition  region  and 
extends  into  the  upstream  solar  wind  and  downstream 
magnetosheath . 

Figure  2  This  figure  depicts  the  planar  bow  shock  model  when 

the  upstream  0gn  is  equal  to  45°.  The  shock 
transition  region  extends  for  100  km  in  the  z 
direction  and  separates  the  upstream  solar  wind  from 
the  downstream  magnetosheath.  The  magnetic  field 
magnitude  and  electron  concentration  are  assumed  to 
increase  linearly  with  increasing  z.  Within  the 


shock  transition  region  the  plasma  flow  velocity  is 

A 

taken  to  be  in  the  z  direction  and  the  plasma  flux  is 
constant.  Due  to  the  planar  geometry,  the  normal 
component  of  the  magnetic  field  is  constant 
throughout  the  transition  region.  A  magnetic  field 
line  is  shown  as  is  6gn  near  the  center  of  the 
transition  region. 

Figure  3  This  figure  depicts  contour  plots  of  the  logarithm  of 

the  model  electron  velocity  distribution  as  a 
function  of  shock  penetration,  z.  The  distribution 
evolves  from  a  relatively  cool  Maxwellian  at  the 
solar  wind  (z  =  0  km)  to  a  relatively  hot  modified 
Lorentzian  at  the  magnetosheath  (z  «  100  km).  As  the 
beam  energy  increases,  the  anisotropy  ratio,  Tj/T|, 
increases  as  does  the  modified  Lorentzian  component 
of  the  model  distribution. 

Figure  4  This  figure  summarizes  the  unstable  frequencies  as  a 

function  of  shock  penetration  for  the  model  distribu¬ 
tion  and  generation  via  the  normal  cyclotron  reso¬ 
nance.  It  is  usually  true  that  for  a  given  frequency 
a  wide  range  of  wave  normal  angles  is  unstable.  The 
increase  in  frequency  with  increasing  z  is  due  to  the 
fact  that  the  normal  cyclotron  resonance  velocity 
must  remain  in  regions  of  large  anisotropy  in  F(z). 


Figure  5  This  figure  summarizes  the  unstable  frequencies  as  a 

function  of  shock  penetration  for  the  model  distribu¬ 
tion  and  generation  via  the  Landau  and  anomalous 
cyclotron  resonances.  It  is  usually  true  that  for  a 
given  frequency  a  wide  range  of  wave  normal  angles  is 

unstable.  The  distribution  is  stable  for  z  greater 

k  u  gp 

than  about  50  km  because  *n — r  — —  evaluated  at  the 

TV 

Landau  resonance  velocity  is  not  sufficiently 
positive. 

Figure  6  This  figure  illustrates  the  relative  orientations  of 

the  magnetic  field  vector,  the  wave  vector,  the  group 
velocity ,  and  the  plasma  flow  velocity  for  generation 
via  the  normal  cyclotron  resonance.  All  vectors  lie 
in  the  plane  of  the  figure.  For  a  given  wave  normal 
angle,  the  two  wave  vectors  that  lie  in  the  plane  of 
the  planar  bow  shock  model  are  distinguished  by  the 
angle  <f>  with  the  <j>  =  0°  direction  defined  to  be 
toward  the  solar  wind.  The  <j>  =  180°  orientation 
lies  on  the  unstable  cone  of  wave  vectors  and  is 
180°  in  azimuth  from  the  <j>  ■  0°  orientation.  For  a 
given  wave  vector  orientation,  the  ray  path  direction 
in  the  shock  frame  of  reference  is  in  the  direction 
of  the  vector  Vg  +  Vf.  For  ky  >  0,  most  of  the  ray 
paths  are  directed  toward  the  solar  wind.  However, 


the  <j>  -  180°  orientation  of  k  together  with  the 
plasma  convection  effect  can  result  in  a  ray  path 
directed  nearly  transverse  to  the  shock  normal. 


Figure  7  This  figure  is  similar  to  Figure  6  but  illustrates 

the  geometry  for  the  case  of  generation  via  the 
Landau  and  anomalous  cyclotron  resonances.  The  ray 
path  directions  in  the  shock  frame  are  given  by  the 
direction  of  the  vector  Vg  +  vf.  In  all  cases,  the 
ray  paths  are  directed  toward  the  magnetosheath. 


Figure  8 
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This  figure  illustrates  the  results  of  the  propaga¬ 
tion  and  integrated  growth  calculations  when  the 
upstream  0Bn  Is  equal  to  45°.  Shown  are  a  magnetic 
field  line  and  several  ray  paths  in  the  shock  transi¬ 
tion  region.  The  ray  paths  for  generation  via  the 
Landau  and  anomalous  cyclotron  resonances  originate 
in  the  lower  left  corner  of  the  figure  while  those 
for  generation  via  the  normal  cyclotron  resonance 
originate  on  the  right  side  of  the  figure.  The  ray 
paths  are  labeled  in  the  format  (up/Zir,  60,  <j>)  where 
0O  and  $  specify  the  Initial  wave  vector  orientation. 
A  $  value  for  0O  -  0°  is  not  applicable  and  so  is 
omitted  from  the  ray  path  labels.  The  solid  portions 
of  the  ray  paths  correspond  to  regions  of  growth 
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Figure  9 
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while  the  dashed  portions  correspond  to  regions  of 
damping.  At  the  heads  of  the  arrows  denoting  the  ray 
path,  the  maximum  integrated  growth  along  the  ray 
path  is  given.  The  ray  paths  terminate  when  the 
integrated  growth  rate  becomes  negative.  It  is  evi¬ 
dent  that  some  of  the  waves  generated  via  the  normal 
cyclotron  attain  large  amplitudes  within  the  shock. 

This  figure  presents  the  results  of  the  propagation 
and  integrated  growth  calculations  when  the  upstream 
0gn  is  equal  to  65°.  The  format  of  the  figure  is 
similar  to  Figure  8.  A  comparison  of  Figure  8  and 
this  figure  suggests  that  the  number  of  wave  packets 
that  attain  large  amplitudes  within  the  shock 
increases  with  increasing  0jjn.  The  largest 
integrated  growth  occurs  for  waves  that  travel 
transverse  to  the  shock  normal. 
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